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1 Introduction 



In 1995 the KARMEN collaboration discovered an anomaly in the time distribution of 
the charged and neutral current events induced by neutrinos from 7r+ and /i"*" decays 
at rest This anomaly may be explained by the production of a new neutral particle 
in pion decay 

7r+^/i+ + x°, (1) 

with the mass 33.9 MeV, barely permitted by the phase space, so that this particle 
moves with non-relativistic velocity. Its subsequent neutrino-producing decays could 
be the source of the delayed neutrinos observed in the experiment. The anomaly 
was recently confirmed by the same group with better statistics and substantially 
reduced cosmic-ray background. 

A search for particles produced in the rare pion decay (|l]) was performed in 
1995 by an experiment at PSI It gave an upper limit for the branching ratio of 
BR(7r+ ^ /i+ + x°) < 2.6 x 10"^ at 95% CL. 

Several candidates for have been proposed in the literature. In ref. 0] the 
authors considered a sterile neutrino, x° = Ug. Their conclusion was that the sterile- 
neutrino hypothesis is compatible with all laboratory constraints, but they noted 
possible problems with astrophysics and cosmology. Further laboratory constraints 
on this model were investigated in ref. ||^ , concluding that mixings of Ug with z/g and 
Ufj^ must be very small, while a mixing with Ur was permitted. In this case Ug would 
predominantly decay through neutral-current interactions into Ur + i + £, where i is 
any light lepton, £ = z/g, z/^, i^r, or e~. The lifetime of Ug with respect to this decay 
was estimated to be in the range 

10"^ sec < r^^ < 150 sec . (2) 

The lower bound on t^^, which comes from the experimental bound on the u^- mass 
(large mixing with Ug makes it too heavy), was discussed in ref. 0. 
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In ref. [0 it was suggested that x° = x could be the hghtest supersymmetric 
particle, photino or zino, that decayed through the channel 

X^l + ^i. ■ (3) 

The new data however, do not agree with the predictions of the model so that 
recently a new version of supersymmetric model was considered [^, according to 
which the light neutralino decayed through a three body channel 

X ^ e+ + e" + z/^,^ . (4) 

In the paper ||^ a new decay mode of muons was proposed as a source of the 
anomaly: 

/i+ e+ + ^ , (5) 
where S* is a scalar boson with mass 103.9 MeV. A search for these decays was reported 



in ref. ||T0| where the upper limit BR(/i'^ ^ e+ + S) < 5.7 x 10 ^ was obtained. This 



limit, though it does not exclude the model, makes it more complicated. 



Recent data by the NOMAD collaboration |TT| permit to strengthen the bound 
on the mixing of Vg with Vr- Their expected lower limit on r^^ is around 0.1 sec. 

On the other hand, cosmology and astrophysics permit to obtain an upper bound 
on Ty^ that may be complementary to direct experiments. If this happens to be the 
case, then the explanation of the KARMEN anomaly by a 33.9 MeV sterile neutrino 



would be ruled out. Recently two papers have appeared WA , where the observation 
of SN 1987A were used to put bounds on the properties of the proposed sterile neutrino 
or light neutralino. In the present paper we derive constraints on the lifetime of a 
33.9 MeV sterile neutrino from both more detailed consideration of SN 1987A and 
of primordial nucleosynthesis. With the present-day accuracy of the data on the 
primordial light-element abundances, the bounds that are found from SN 1987A tend 
to be stronger than those found from primordial nucleosynthesis. However, the latter 
remain interesting, first, because they may be competitive in the nearest future with 



an improved accuracy of the data on primordial nucleosynthesis, and second, in a 
hypothetical case where Ug would have an anomalously strong (stronger than the 
usual weak) interaction with nucleons. In that case the SN 1987A bounds would not 
apply, while the nucleosynthesis bounds would survive. 

2 Primordial Nucleosynthesis 
2.1 General Features 

A heavy unstable sterile neutrino would influence big-bang nucleosynthesis (BBN) 
through its contribution to the cosmological energy density by speeding up the ex- 
pansion and enlarging the frozen neutron-to-proton ratio, r„ = n/p, and less directly, 
though stronger, through its decay products, z/g, z/^, and Ur- The impact of z/^ and Ur 
on BBN is rather straightforward: their energy density increases with respect to the 
standard case and this also results in an increase of r„. This effect can be described 
by the increased number of effective neutrino species N^, during BBN (in the standard 
case A^",^ = 3). The increase of the energy density of z/g, due to decay of z/^ into z/g, 
has an opposite effect on r„. Though a larger energy density results in faster cooling, 
the increased number of Uf. would preserve thermal equilibrium between neutrons and 
protons for a longer time and correspondingly the frozen n/p-ratio would become 
smaller. The second effect is stronger, so the net result is a smaller n/p-ratio. There 
is, however, another effect which is related to the distortion of the energy spectrum 
of z/e from the decays of z/^. If the spectrum is distorted at the high-energy tail, as is 
the case, then proton formation in the reaction n + Uf. p + e~ would be less efficient 
than neutron creation in the reaction + p ^ n + . We found that this effect is 
quite significant. The decays of z/^ into the e"''e~-channel will inject more energy into 
the electromagnetic part of the primeval plasma and this will diminish the relative 
contribution of the energy density of light neutrinos and diminish r„. 

We could use the technique and numerical code from our earlier papers where the 
effects of non-equilibrium massless neutrinos (in the standard model) fl^ and possibly 
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massive [|T5] and/or unstable |jT6| tau- neutrinos were precisely calculated. However, 



since we do not need the accuracy of a fraction of per cent achieved in these papers, 
we will use a simpler and considerably less time consuming approximate approach 
that we will now discuss. 

2.2 The model 

If the KARMEN anomaly is explained by a heavy sterile neutrino with m^^ = 
33.9 MeV then, as was mentioned in the Introduction, it may only mix with u^-, 

Vr = T-'i cos 9 + z/2 sin 9 , 

Vs = —J^i sin 9 + z/2 cos 9 , (6) 

where 9 is the vacuum mixing angle and z/i and z/2 are the mass eigenstates; the mass 
difference is positive: Sm^ = m\ — m\ (33.9 MeV)^ > 0. Through neutral-current 
interactions Vg could decay into Vj. and a pair of other light leptons. The corresponding 
processes and their matrix elements are presented in table 1. The lifetime of is 
given by the expression 
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5.7 X 10-Sec 
(sin2 29)/4 



(7) 



where 

= — 1/2 + sin^ ^vy and (yfR = sin^ ^VK • (8) 
According to the combined experimental data, the lifetime lies in the range 

0.1 sec < Ty^ < 150 sec (9) 

and correspondingly 

3.9 X 10"^ < sin 29 < 0.15 . (10) 

Even with a very small mixing, z/^ could be abundantly produced in the early universe 
when the temperature was higher than niy^. Their production rate can be estimated 



as Bin 



r,, = ^ sin^ 2eMrH., (ii) 

where F^y = 2.5 G^pT^ is the averaged weak interaction rate and Gm is the mixing 

angle in the medium. According to the calculations of refs. ||T7| , |19| 

sin 2eM ~ isin2e)/2 ^ (^in2e)/2 

^ 1 + 0.76 X 10-19 T6(5m2)-i 1 + 6.6 X 10-23 T6 ' ^ ^ 

where T and Sm^ are taken in MeV. One sees that matter effects are not important 

for T < 5 GeV. 

Comparing the production rate ([TT|) with the Hubble expansion rate 



we find that sterile neutrinos could be abundantly produced in the early universe if, 
roughly speaking, (sin20)2(T/3 MeV)^ > 1. Even for G at the lower limit given by 
the relation (p!0D, the equilibrium condition is fulfilled for T ^ 120 MeV. For this 
small value of G the original equilibrium number density of Ug would be diluted at 
smaller T by annihilation of pions and muons. However, as we see in what follows, nu- 
cleosynthesis strongly disfavors large values of t^^. Correspondingly, for large mixing 
angles Ug remains in equilibrium for much smaller T, and this dilution is not essential. 

The evolution of Us at lower temperatures, T < rriiy^, is considered below. We 
calculate the number density of the heavy Us assuming that initially they were in 
equilibrium and in the process of freeze-out they interacted with the thermal equi- 
librium bath of light particles. With the evident correction for the Ug decay, these 
calculations are very similar to the usual freeze-out calculation of massive species. 
There is one important difference however. Normally massive particles disappear in 
the process of mutual annihilation, so that the rate of freezing is proportional to the 
number density of the particle in question, hm/rim ~ cra,nnnm, which becomes expo- 
nentially small when T < m. Sterile neutrinos may also disappear in collisions with 
massless leptons 

i^s + h-^h + h, (14) 



so their extinction through this process would be more efficient than by mutual anni- 
hilation, Us + i^s —>■ all- The decoupling temperature of the process ([I^) can be approx- 
imately found from the decoupling of the usual massless neutrinos at ~ 2 MeV 
by rescaling it by the mixing parameter, (sin2B)^/'^. To obtain a better estimate 



we solve the corresponding kinetic equation for Ug in sec. p.3| . This permits us to 
determine the distribution function, f^^{x,y), where 

X = rrioait) and y=pa{t), (15) 

are convenient dimensionless variables with a{t) the cosmic scale factor and tuq the 
normalization mass that we choose as mo = 1 MeV. In what follows we will often 
skip mo, keeping in mind that the relevant quantities are measured in MeV. In terms 
of these variables, the kinetic equations have the form 

{dt - Hpd,)f = HxdJ = /,oii , (16) 



where 



'coll 
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2E J Y \'2Ei{2Tr)^ J y V2^/(27r) 
X (27r)V^)(Ep.-EP/)l^^/l'W.//) , (17) 

i f 

is the collision integral and 

m, //) = - n n(i - //) + n // n(i - /^) (is) 

if f i 

with sub-i and sub-/ meaning initial and final particles. 

After the distribution function of is found, the next step is to find the distri- 
bution functions of the light neutrinos and in particular their energy densities. The 
distributions of electrons and positrons are of course assumed to be very close to 
equilibrium because of their very fast thermalization due to the interaction with the 
photon bath. However, the evolution of the photon temperature, due to decay and 



7 



annihilation of the massive Ug, becomes different from the standard one, ~ 
by an extra factor (1 + A) > 1: 

= [1 + A{x)]/x . (19) 

At sufficiently high temperatures, T > ~ 2 MeV, light neutrinos and elec- 
trons/positrons were in strong contact, so that the neutrino distributions were also 
very close to the equilibrium ones. If disappeared sufficiently early, while thermal 
equilibrium between and neutrinos remained, then i/g would not have any ob- 
servable effect on primordial abundances, because only the contribution of neutrino 
energy density relative to the energy density of and 7 is essential for nucleosyn- 
thesis. Hence a very short-lived i^g has a negligible impact on primordial abundances, 
while with an increasing lifetime the effect becomes stronger. Indeed at T < the 
exchange of energy between neutrinos and electrons becomes very weak and the en- 
ergy injected into the neutrino component is not immediately redistributed between 
all the particles. The branching ratio of the decay of i/g into e'^e~ is approximately 
1/9, so that the neutrino component is heated much more than the electromagnetic 
one. As we mentioned above, this leads to a faster cooling and to a larger n/p-ratio. 

So, for the numerical calculations we adopt the following procedure. We assume 
that at sufficiently high temperature, e.g. — b MeV, or equivalently at the initial 
value of the scale factor Xi — 1/Tj = 0.2, there is a complete thermal equihbrium 
between active neutrinos and electrons. Starting from that moment we calculate the 
corrections to the active neutrino distribution functions 

5U = U- , (20) 

where the equilibrium distribution function is assumed to have the standard Fermi- 
Dirac form with a temperature that drops as 1/x, 

/,^^ = (e^ + l)-^ . (21) 

The evolution of the photon temperature, i.e. A(a;), is determined from the energy 
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balance equation 

dp/dx = -•i{p + p)/x , (22) 

where p and p are respectively the total energy and pressure densities in the cosmo- 
logical plasma. 

It is worth noting that we normalized the scale factor in such a way that at the 
initial moment XjTj = 1. If we change the initial moment for calculating A, it would 
result in a different definition of x, but this is unobservable. We have checked that for 
Xi = 0.1-0.3 the results weakly depend upon the value of Xi, and that the corrected 
neutrino distribution fi, = f^^ + 5f quite accurately maintain the equilibrium shape 
with the same temperature as electrons, = at these early times. 

We made several assumptions that permitted to simplify calculations very much: 
Us was assumed non-relativistic, the equilibrium distribution functions "inside" the 
collision integral were taken in the Boltzmann approximation, while "outside" they 
were taken in Fermi-Dirac form, and the kinetic equations for 6fi, were taken in a 



simplified form, so that some fraction of neutrino energy was lost (see sec. p.4|) . All 
these assumptions lead to a weaker impact of Us on nucleosynthesis, so the real bound 
should be somewhat stronger than what is presented below. 

One more comment is in order. We do not take into account oscillations between 
Us and Ur for T < m^^. This is perfectly justified because in the interesting range of 
neutrino energies the oscillation frequency is so high and the velocities of z/^ and Uj. 
are so much different that the coherence is quickly lost and they can be considered 
as independent particles. Medium effects are also not important for the considered 
positive mass difference ^10^ MeV^. 

2.3 Evolution of Heavy Neutrinos 

The evolution of the occupation number of z/^ is determined by its decays and inverse 
decays, listed in table 1, and by the reactions (|1^ with all possible sets of light leptons 
permitted by quantum numbers, presented in table 2. Taking both contributions into 
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the collision integral and assuming that the massless species are in thermal equilibrium 
with temperature T and that both helicity states of Ug are equally populated, we 
obtain 

1.48 X / 10.75 V^' /.I"-/.. 



+ ^f^ + TTT + 



E,^ ' ml \ 4 144 V ml SE.^ml 



(23) 



where Ei,^ = ^Jm^ + (y/x)'^ and p^^ = y/x are the energy and momentum of Ug 
respectively and g*(T) is the effective number of massless species in the plasma deter- 
mined as the ratio of the total energy density to the equilibrium energy density of one 
bosonic species with temperature T, g^: = ptot/(7r^T'^/30). The coefficient 1.48 comes 
from expressing the lifetime in sec according to the relation 1/MeV = 0.658 x 10"^^ sec 



and from the value of the Hubble parameter H = y 87rptot/3M|[. The first term in 
square brackets, m^jE,^^, comes from summing the squared matrix elements in ta- 
ble 1 (decay) while the rest is obtained by using that the sum of the squared matrix 
elements from table 2 (collisions) can be written as 

|M|2 ^{l+gl + gl) [(pi . p2)(p3 " Pa) + 2(pi ■ p,){p^ ■ P2)] • (24) 

The assumption of Boltzmann statistics for massless leptons considerably simpli- 
fies the calculations, because the integral over their momenta can be taken explicitly. 
This approximation results in a larger value of the collision integral, i.e. in faster 
decay and reaction rates and correspondingly to a smaller abundance of Vg. Thus the 
restriction on r^^ obtained in this approximation is weaker than the real one. 

At high temperatures, T > 5 MeV, (^) was integrated with the simplifying 
assumption Tx = 1. This also leads to a weaker bound on t^^. Assuming entropy 
conservation one can check that Tx may change by a factor 1.1 due to Ug decays and 
annihilation. In fact the effect is stronger because the number density of Ug is much 
larger than the equilibrium one. However, since we started at T = 5 MeV < m^^, 
when the energy density of Ui,^ was already somewhat suppressed, the variation of Tx 
from that moment would be weaker. 
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Naively one would expect a variation of f^^ similar to the variation of Tx. However, 
it is argued below that the effect may be much stronger. Indeed, since z/^ disappears 
in the collisions with massless particles, its number density is much more sensitive 
to the variation of the coefficient in front of the collision integral. In the case of the 
standard freezing in two-body annihilation the frozen number density, nj, is known 
to be inversely proportional to the annihilation cross-section. This is not true for the 
case of reactions (0). One can solve the kinetic equation explicitly and check that 
the result is exponentially sensitive to the coefficient in front of (/^^ — fuj- Therefore 
we took the variation of Tx into account, starting from Xj = 0.2, according to the 
energy balance law (p2D, see also (p4D. A simpler and more common method based 
on entropy conservation is not accurate enough because the sterile neutrinos strongly 
deviate from equilibrium in the essential range of temperatures and thus entropy is 
not conserved. 

In the non-relativistic limit we may integrate both sides of (^) over d^y /{27c)^ to 
obtain the number density of Ug in the comoving volume 



n,. 



3 



where 

B{x) 



lA8x 



10.75 



ctxi— exp 
dxi 



PX 

/ dx2B{x2) 

J XI 



1/2 



(Tx)2(r,7sec) \g,{T) ^ 



3x2'{Txr /3C(3) , ^(Tx 



x^mf. 



+ 



xrrii. 



(25) 



(26) 



Usually Tx is assumed to be constant, quite often normalized as Tx = 1. A small 
variation of this quantity usually is not very important for the value of n. In our 
case, however, estimating the integral by the saddle point method reveals that the 
dependence on Tx appears in the exponent and the effect may be significant even for 
a small variation of Tx. The dependence on Tx for the scattering term results in an 
increase of fij^^, while that for the decay term results in a decrease of Uiy^. 

The results of the numerical solution of (p3D for different values of lifetimes are 
presented in figs. ^ ^ and |^. The first one shows the evolution of x^riiy^ for t,^^ = 
0.1, 0.2, and 0.3, the second is x^p^^, while the third figure is a snap-shot of the 
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distribution functions, f„^, at x = 1. All the distributions are significantly higher 
than the equilibrium ones. For example, the equilibrium distribution is about 10 
orders of magnitude below the calculated curve for t^^ — 0.1. 

2.4 Light Neutrinos 

At the onset of nucleosynthesis, practically disappeared from the primeval plasma, 
at least for sufficiently small lifetimes, which are near the bound obtained below. 
However, the products of their decays and annihilation distorted the standard nucle- 
osynthesis conditions: the energy density as well as the spectrum of neutrinos were 
different from the standard ones, and this changed the light element abundances. 

We will look for the solution of the kinetic equations governing the evolution of 
the distribution functions of light neutrinos in the form 



where a — e, or r. The first term is equal to the equilibrium distribution function 
for the case when the temperature drops as the inverse scale factor, i.e. Tx = 1. In 
reality this is not the case and the function 6f takes into account both the variation 
of Tx and the spectrum modification of the active neutrinos. 

The distribution of always has the equilibrium form, f^ — [l + exp(£'/T)]~^ = 
[1 + exp{y/Tx)]''^ and the product Tx is taken in the form 



where A is assumed to be small, such that a perturbative expansion in A can be 
made. A similar first-order perturbative expansion is made with respect to Sf, so 
that the collision integral becomes linear in terms of A and Sf. 

It is assumed that initially both 5f and A were zero. To this end one should find 
an appropriate value of the initial "time", Xi, or temperature, Tj. The temperature 
should be sufficiently small, such that the number density of is already low and 
hence the function A would not rise too much. On the other hand, the temperature 



(27) 



= 1 + A{x) , 



(28) 
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should not be too low, otherwise equilibrium between neutrinos and electrons would 
not be established and the assumption of 5/ = could be grossly wrong. As we see 
in what follows, convenient initial values that satisfy both conditions are Xi ~ 0.2 
and Tj ~ 5 MeV for t^^ ~ 0.2 sec. Of course, the later rise of temperature by the well 
known factor 1.4 because of e+e~-annihilation is taken into account explicitly. 

2.5 Collision Integral and Source Term 

There are two kinds of terms in the collision integral, "easy" ones where the unknown 
function 6f depends upon the external variable y, which is not integrated upon, and 
"difficult" ones when 6f is under the integral. The terms of the first kind come, with 
negative sign, from external particles in the initial state. They give 

/10 7'^\^^^ 

dJf.Ax, y) = -0.26 i-j-j 5/.,(l +gl + gl){y/x') , (29) 
where gi = sin^ 9w + 0.5, while gji = sin^ 9w For z/^ and the functions g^ and gji 



are given by expressions (^) . Using the simple expression ( PD| ) one finds that the rate 
of approach of Ue to equilibrium is given by 

5/,^ ~exp(-0.13y/a;3) . (30) 

One sees that equilibrium would be efficiently restored for i/g if x < 0.5 y^^^ or T > 
2?/^^/^ MeV, while for u^^j- equilibrium would be restored if T > 2.25 y^^^^ MeV. 
These results are close to the standard estimates known in the literature (for the 



most recent ones see e.g. 1^). However, there are the following drawbacks in this 
derivation. First, it was done under the assumption that the source of distortion acted 
for a finite time and thus (|30|) could be valid only after the source has been switched 
off. If there is a constantly working source, S{x,y), equilibrium is always distorted 
roughly by the factor S/T, where F is the effective reaction rate. Another and more 
serious argument against the validity of (|30| ) is that the "difficult" part of the collision 
integral was neglected. One can see from the general expressions ( [T7|JT^ ) that 6f also 
appears under the integral over momenta, and it can be seen that it comes mostly with 
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a positive sign. These terms counteract the smoothing action of expression (p9|) and 
shift equihbrium restoration to considerably higher temperatures. One can take them 
into account exactly, numerically solving the kinetic equations with the exact collision 
integral that can be reduced down to a 2-dimensional integration over particle energies 
(see e.g. ^^), but this is a very time consuming procedure. Instead we will here use 
a much simpler approach. We will approximately represent such "difficult" terms by 
the integrals over energy chosen in such a way that the kinetic equations satisfy the 
law of particle conservation in the comoving volume if only elastic scattering is taken 
into account, i.e. the kinetic equation should automatically give dx{x^n) = 0. The 
exact equations should simultaneously satisfy energy conservation law, but working 
with the exact equations is much more complicated. We checked the validity of our 
approximate procedure for a case of a more simple reaction amplitude, where we 
compared the approximate and exact results and found very good agreement. The 
approximation that we use breaks the energy conservation law in the kinetic equation, 
so that some small fraction of energy going into light particles from Ug decays and/or 
annihilation is lost. This diminishes the effects that we are discussing so that the real 
constraints on t^^ should be stronger. 

The kinetic equations for light neutrinos in this approximation can be written as 

1/2 

9* 



dJUAx.y) = 5..(x,i/) + 0.26 f {l + gl + gl){y/x') (31) 



x{-^/. + ^^^^[l + 0.75(.i + .^)] 



dy2ylSU^{x,y2) + I I dy2yl {6U^{x,y2) + 6f^^{x,y2) 



15 1 + gt + gj, 

X di 

The coefficients (^i^/j for i/^,- are given by (|D and for v^, are presented after (|29|). 
The source term S describes injection of non-equilibrium neutrinos by z/g decays or 
reactions (0) with light leptons. In what follows we include only decays. We have 
estimated the contributions of the reactions (0) and found that they slightly improve 
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the restrictions that we have obtained. For z/g, t'^, and Vr the contributions of the 
decay term are respectively 



0.024 / 10.75 



1/2 



+ 77(1 + ^!+^?^ ! 



X (n,^-<^)^(m,,x/2-y) , (33) 



where is the number density of Vg and 9{y) is the step function which ensures 

energy conservation in the decay. The factor 0.012 comes from the product of the 
branching ratio BR = 96/4/(1 +gl + gji) = 21.3 with the factors 3/m^^ = 7.7 x 10^^ 
{rriiy^ in MeV) from the normahzation, we divide by = 0.203 from the number 
density, and 6.582 from the relation between MeV and sec. Dividing by the Hubble 



parameter gives a factor, 1. 221/(1. 88^87r/3) = 0.2244, and we find 21.3 x (7.7 x 
10-^) X 0.2244 X 6.582/0.203 = 0.012. 

The coefficient in front of the collision integral is in reality momentum dependent, 
and hence is slightly different from 0.4 for (and 0.29 for t'^,r), which are often 
used in the literature. We extracted the correct momentum dependent coefficients 



from our Standard Model code ||Tj], and used this in the calculations. This slightly 
weakens the bound on the lifetime. 

The function A(x) is determined from the energy balance condition (|2^) which in 
the present case reads 



dA 



dx Ax^p 



EM 



xd{x^p^J ^ d{x^5py) 



dx dx 



(34) 



where pem is the energy density in the electromagnetic sector, and p^^ and p^ are the 
energy density of Vg and the total energy density of all light neutrinos respectively. We 
also used another method for calculating Tx starting from earlier times and obtained 
stronger results, so we believe that the limits that we obtain here are quite safe. 
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2.6 A Few Numerical Technicalities 

Let us summarize a few technicalities related to the numerical approach. We divide 
the time into 3 regions. First we integrate only (^) from very high temperatures, 
T = 50 MeV {x = 0.02) and until T = 5 MeV (x = 0.2). We assume that initially the 
sterile particles are in equilibrium (see sec. p^.2| ). In this way we can follow the freeze 
out and initial decay of the sterile particles, without worrying about the equilibrium 
active neutrinos and the electromagnetic plasma. In the next region, 0.2 < x < 50, 
we solve (|23| , |34D with the initial values A = and 5fv^^^^^ = 0. It is worth 



mentioning that from the very beginning we separately calculate and include the 
annihilation of the electrons, which increase the photon temperature with a standard 
factor 1.4. Finally, for very high x > 50 (T < 0.03 MeV) when all sterile neutrinos 
have disappeared and the active neutrinos have long decoupled, we solved only the 
kinetic equations governing the ra-p-reactions needed for the nucleosynthesis code. We 
use an 800 point grid in momentum in the region < y < 80, and checked that the 
results are insensitive to the doubling of the grid. For the BBN calculations we use 
rjw = 5. 

2.7 Results 



We have solved (|23| , |34|) numerically for different lifetimes. In fig. |l| we plot the 
evolution of the normalized number density, function of x. One sees for 

T^^ = 0.3 how Us first freeze out, followed by the subsequent decay. In fig. || a similar 
plot of the normalized energy density, x^p^^, is presented as a function of x, and the 
effect of the sterile particle being massive is evident. In fig. ^ we present a snap-shot 
of the distribution function at the time x = 1. The equilibrium distribution function 
at this time is about 10 orders of magnitude smaller than the curve for r = 0.1 sec. 

The calculated energy densities of all light neutrino species relative to the elec- 
tromagnetic energy density, pe + P75 are presented in fig. ^ as functions of x. One 
sees that this fraction is higher for longer lifetimes, especially around the time x = 1, 
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when the n/p-ratio freezes out, leading to an expected increase in the final helium 
abundance. In fig. § a snap-shot of the spectrum of z/g, namely y'^fu^ and the distor- 
tion y^Sfu^ are presented for x = 1 for the lifetimes 0.1,0.2,0.3. A distortion of the 
electronic neutrino spectrum has a strong impact on nucleosynthesis, while and 
Ut- act only by their total energy density. It is noteworthy that the increase in p^y^ 
acts in the opposite direction to an increase in Pu^^vri since it reduces the effective 
number of light neutrinos, or in other words, it gives rise to a smaller mass fraction 
of primordial ^He, while an increase in the high energy part of z/g spectrum results in 
a larger mass fraction of ^He. 

The results of the calculations have been imported into the modified Kawano |pT| 
nucleosynthesis code, and the abundances of all light elements have been calculated. 
At each time step x, we find the corresponding photon temperature and total energy 
density. Furthermore we integrate the kinetic equation governing the n/p evolution 
taking into account the distorted spectrum of v^- The final helium abundance is 
presented as a function of the z/^ lifetime in fig. ^ By translating these results into 
effective number of neutrinos one sees that if we allow for AA^ = 0.2 (as suggested 
by p2|), then only lifetimes lower than t^^ = 0.17 sec are permitted. If one is more 



conservative and allows for one extra neutrino species, AA^ = 1.0, then lifetimes 
longer than r^^ = 0.24 sec are excluded. A drop in the helium abundance around 
Tj,^ = 0.1 is related to the dominant role of the energy density, since the spectrum 
distortion is shifted to smaller energies. 

Finally in fig. |^ we compare the KARMEN experimental data with the bound 
obtained from BBN. One sees that the expected new data from the NOMAD Col- 
laboration together with our BBN bound leave a small allowed window for a sterile 
neutrino with a lifetime around 0.1-0.2 sec. 

3 Supernova Limits 
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3.1 Small Mixing Angle 

The mixing angles of sterile neutrinos with the standard active flavors is tightly 
constrained by standard arguments related to supernova (SN) physics and to the 
neutrino observations of SN 1987A. Some of these arguments have been sketched out 
in the context of the KARMEN anomaly in Ref . @] . 

The simplest limit arises from the "energy-loss argument." The SN 1987A obser- 
vations imply that a SN core may not emit too much energy in an "invisible channel" 
as this would unduly shorten the observed neutrino burst. Reasonably accurate limits 
are obtained by demanding that the "exotic" energy- loss rate should obey p3| , p4 | 



e 



< 10^^ erg g-i , (35) 



where e is to be calculated at typical average conditions of a SN core (p = 3 x 
10^^ g cm-3, T = 30 MeV). 

Sterile neutrinos Ug are produced because they mix with one of the standard ones. 
If the standard neutrino masses are all in the sub-eV range, the assumed sterile 
neutrino mass of m^^ = 33.9 MeV assures a mass difference so large that medium 
effects on the oscillations can be neglected. The oscillation frequency is so large that a 
standard neutrino Ua, once produced, oscillates many times before collisions interrupt 
the coherent development of the flavor amplitude. The average probability of flnding 
the original z/^ in the z/^ flavor state is |sin^(2G) where B is the z/a-i^g-mixing angle. 
Oscillations are interrupted with the collision rate F of the z/^ flavor, leading to the 
standard sterile-neutrino production rate of | sin^(26)r. 

We estimate F as the neutral-current collision rate on free nucleons, ignoring 
correlation and degeneracy effects, so that: 

„ CI + 3C^ 



GIubEI , (36) 



TT 



with Gp the Fermi constant and ns the number density of baryons (nucleons). For a 
mix of protons and neutrons we use an average neutral-current coupling constant of 
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{Cy + 3C^) ~ 1. If we further assume that the trapped active neutrinos do not have 
a significant chemical potential (true for z/^ and Ur, but not for z/g), the energy- loss 
rate is 

sin|e)^,.«^^^^ (37) 

2 TT'^mj^f Jo e^-"'^ + 1 

where is the nucleon mass. For simplicity we use Maxwell-Boltzmann statistics 
for the neutrinos (we ignore the +1 in the denominator under the integral) and find: 

sin2(2G) 120 G^T^ 
= sin2(2e) 2.8 x 10^^ erg g'^ s^^ T^^ 

(38) 

where T30 = T/30 MeV. Comparing this result with (|35|) leads to a limit 

sin2(2e) < 3 X 10"^ Tg-^ . (39) 

The temperature T = 30 MeV is at the lower end of what is found in typical numerical 
calculations so that this limit is reasonably conservative. 

For the mixing with Uf, the limit is more restrictive because the electron neutrinos 
in a SN core are highly degenerate, leading to a larger conversion rate and a larger 
amount of energy liberated per collision. Several authors found sin^(2B) ^ 10^^° for 
this case p5|, ^ . 



In our calculation of the emission rate we have taken the Ug to be effectively 
massless. The average energy of trapped standard neutrinos is about 3T which far 
exceeds m^^. The scattering cross section scales with E"^, favoring the emission of 
high-energy neutrinos. The average energy emerging from a SN core is thus found 
to be about 5T so that neglecting m,^^ is a good approximation — the Ug are highly 
relativistic. 

The KARMEN experiment implies that the lifetime of Ug exceeds about 0.1 s so 
that these particles escape from the SN before decaying. On the other hand, they 
must decay on the long way from SN 1987A to us. If the decay products include 
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e'^e pairs one will also get 7 rays from inner bremsstrahlung in the decay. The 
non-observation of a 7 ray burst coinciding with SN 1987A then leads to further 
limits ^ g. 

The decay products likely would include standard neutrinos which would have 
shown up in the detectors. However, because of the high energies of the sterile 
neutrinos which are representative of the SN core temperature, such events would 
have much larger energies than those expected from thermal neutrino emission at the 
neutrino sphere. Therefore, the emission of sterile neutrinos from the core and their 
subsequent decay cannot mimic the standard SN neutrino signal. 

3.2 Large Mixing Angle 

These upper limits on the mixing angle are only valid if the sterile neutrinos actually 
escape from the SN core after production, a condition that is not satisfied if one of the 
mixing angles is too large. Since the mixing angle between z/g or z/^ and Vg is already 
constrained from laboratory experiments to be small in this sense, we worry here only 
about the z/^-z/^ mixing angle. We continue to assume that the standard-neutrino mass 
eigenstates are small so that the mass difference and hence the oscillation frequency 
between z/^ and Vg will be large compared to a typical collision rate in a SN core 
and that the mixing angle in the medium is identical to the vacuum mixing angle. 
Therefore, after production the chance of finding a Vg in the active Vj. state will be 
given by the average value |sin^(20). If the mixing angle is large this will mean 
that the sterile neutrino essentially acts as an active one. It will be trapped, and the 
energy- loss argument is not applicable. 

If the mixing angle is not quite maximal, sterile neutrinos will still be trapped, but 
their mean free path will be larger than that of an active flavor. Energy is transported 
out of a SN core by neutrino diffusion, a mechanism which is more effective if the mfp 
is larger. Simply put, a neutrino can transport energy over distances of order the mfp 
so that more distant regions are thermally coupled if the mfp becomes larger. Most of 
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the transporting of energy is done by the most weakly coupled particles which are still 
trapped. For example, in a SN core the photon contribution to the energy transport 
is negligible because their mfp is very much smaller than that of standard neutrinos. 
Conversely, the contribution of sterile neutrinos will increase with an increasing mfp, 
i.e. with a decreasing sin^(20). 

The effect on the SN 1987A signal will be identical to the effect of freely escaping 
sterile neutrinos, i.e. the signal will be shortened. We stress that the signal duration 
is determined by the diffusion time scale throughout the star. Therefore, increasing 
the mfp in the deep interior of the star shortens the cooling time scale. In a numerical 
study the neutrino opacities were artificially decreased. It was found that the 
efficiency of neutrino transfer in the star should not be more than about twice the 
standard value to remain consistent with the SN 1987A signal characteristics. Like- 
wise, in a different study the number of standard neutrino flavors was artificially 



increased, again leading to an increased efficiency of energy transfer. Doubling the 
effective number of neutrino flavors appears excluded from the SN 1987A data. 

The interaction rate of sterile neutrinos is that of a standard Ur, times | sin^(20). 
Conversely, the collision rate for i/^ is 1 — | sin^(2G) times the standard rate because 
a Uj- has an average chance 1 — | sin^(2B) of being measured as a Ug. Assuming that 
the standard transfer of energy is dominated by u^, and their anti-particles (the 
mfp for Vf. and is much shorter due to charged-current reactions), then adding the 
sterile neutrino will enhance the rate of energy transfer by a factor 

ITT ^ 1 - isin2(2e) + isin2(2e)) ■ 

Maximal mixing corresponds to sin^(20) = 1, implying that both Vr and each 
scatter with half the standard rate, so their mfp is each increased by a factor of 2. 
Moreover, the sterile neutrino contributes a second channel for the transfer of energy. 
Therefore, the energy flux carried by maximally mixed Vs and Vr is four times that 
carried by a standard Vr. This explains the limiting behavior of (|40D for maximal 
mixing. 
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The new experimental limit on the z/^-z/^ mixing angle is sin^(2G) < 10~^. There- 
fore, the efficiency of energy transfer out of a SN core would be enhanced by more 
than two orders of magnitude. Such an enhancement is certainly not compatible with 
the SN 1987A signal, implying that the mixing angle has to be very small, i.e. that 
it must obey (|39|). Therefore, the loop hole of a large i^^-z/^ mixing angle has been 
plugged by the new experimental constraints. 

4 Conclusion 

Cosmological and astrophysical arguments seem to exclude the interpretation of the 
KARMEN anomaly by an unstable sterile neutrino mixed with v^-. The arguments 
based on SN 1987A are stronger than the nucleosynthesis bound. The supernova limit 
is given by (P^D, i.e. t^^ > 6 x 10^ sec. This result, together with the direct experimen- 
tal constraints on v^-, completely excludes a 33.9 MeV sterile neutrino. Primordial 
nucleosynthesis permits to exclude roughly > 0.2 sec. So with the existing direct 
experimental limits, some window, t^^ = 0.1 — 0.2 sec, remains open. The simplifica- 
tions that we made in deriving the BBN bound typically lead to a weaker result, so 
the real bound may be somewhat stronger. However, if exact calculations conffim the 
decrease of ^He around t,^^ = 0.1 found in this paper, then BBN will never exclude 
lifetimes ~ 0.1 sec. 

The effect happens to be surprisingly sensitive to usually neglected phenomena, 
particularly the shape of the z/g spectrum and to the non-adiabatic variation of tem- 
perature. We hope to do exact (but rather long) calculations later. Together with 
a possible improvement of the observational data on light element abundances and 
a better understanding of the subsequent changes of light elements in the course of 
cosmological evolution, the BBN bound may become competitive with the supernova 
one. If so, the very attractive hypothesis of a 33.9 MeV sterile Ug would be killed by 
two independent arguments. Moreover, one may imagine a case when the supernova 
arguments are not applicable, while the nucleosynthesis ones still operate. For exam- 
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pie, if Us possesses an anomalously strong interaction with nucleons (stronger than 
the usual weak one), then it would not noticeably change the energetics of supernovae 
that we described here, but would affect nucleosynthesis practically at the same level 
as discussed in section 0. This new interaction might be related to the anomalously 
high mass of z/^. If this is the case then a small window for a 33.9 MeV Ug with 
lifetime r^^ = 0.1-0.2 sec, still exists at the present time, and stronger experimen- 
tal bounds, as well as more accurate calculations of the impact of Ug on primordial 
nucleosynthesis, are needed. 
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Table 1: Matrix elements for decay processes; 9l = 9l — ^ = —\ + siv? 6w and 
qr = sin^ 9w 
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Table 2: Matrix elements for scattering processes; gL — gL — i — —^ + sin^ 9w and 
qr = sin^ 9w 
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Figure 1: The normalized number density, function of x for the hfetimes 

0.1 sec (sohd), 0.2 sec (dashed) and 0.3 sec (dotted). 




Figure 2: The normalized energy density, x"^ p^^, as a function of x for the hfetimes 
0.1 sec (sohd), 0.2 sec (dashed) and 0.3 sec (dotted). 
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Figure 3: A snap-shot of the distribution function, f^^ {y) , at the time a; = 1 for the 
hfetimes 0.1 sec (sohd), 0.2 sec (dashed) and 0.3 sec (dotted). 




Figure 4: The energy density of aU the active neutrinos divided by the energy density 
in the electromagnetic plasma, J2 Pu/ Pem, as a function of x for lifetimes 0.1 sec 
(solid), 0.2 sec (dashed) and 0.3 sec (dotted). 
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Figure 5: A snap-shot of the spectrum of z/g, namely y'^fi,^ and the distortion y^^jv^ 
at a; = 1. The hfetimes are 0.1 sec (sohd), 0.2 sec (dashed) and 0.3 sec (dotted). 
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Figure 6: The final helium abundance as a function of lifetime. The horizontal lines 
correspond to AA?" = 0.2 and 1.0 extra effective neutrino species. 
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Figure 7: Branching ratio versus lifetime. A comparison of the KARMEN experi- 
mental data and the bound obtained here. The BBN bound excludes lifetimes bigger 
than 0.17 sec. The SN 1987A bound excludes the entire region presented. 
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